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Abstract. Spectroscopic observations covering the wavelength range 3600-4600A are presented for a large flare 
on the late type M dwarf AT Mic (dM4.5e). A procedure to estimate the physical parameters of the flaring 
plasma has been used which assumes a simplified slab model of the flare based on a comparison of observed and 
computed Balmer decrements. With this procedure we have determined the electron density, electron temperature, 
optical thickness and temperature of the underlying source for the impulsive and gradual phases of the flare. The 
magnitude and duration of the flare allows us to trace the physical parameters of the response of the lower 
atmosphere. In order to check our derived values we have compared them with other methods. In addition, we 
have also applied our procedure to a stellar and a solar flare for which parameters have been obtained using other 
techniques. 
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1. Introduction 

Stellar flares are events where a large amount of energy is 
released in a short interval of time, radiating at almost 
all frequencies in the electromagnetic spectrum. Flares 
are believed to result from the release of magnetic en¬ 
ergy stored in the corona through re connection (see re¬ 
views by Mirzoyan 1984 Butler 1991 ; Haisch et al. 1991 


Garcia Alvarez 2000|). Many types of cool stars produce 


largest solar flares observed involve energies of ~ 10 32 erg 
( pershberg 1989| ), while large flares on dMe stars can be 
one order of magnitude larger (Doyle & Mathioudakis 
199C; Byrne fc McKay 1990|). Even more energetic flares 


occur on the RS CVn binary systems, the total flare en¬ 
ergy in the largest of this type of systems may exceed 
E~ 10 38 erg (Doyle et al. 1992; Foing et al. 1994). Such 


a change in the star’s radiation field modifies drastically 
the atmospheric properties over large areas, from photo- 


spheric to coronal layers. Models by Houdebine (1992) in¬ 


flares ( Pettersen 1989| ), sometimes at levels several orders 
of magnitude more energetic than their solar counterparts. 
The exact mechanism(s) leading to the energy release and 
subsequent excitation of various emission features remains 
poorly understood. In the dMe stars (or UV Cet type 
stars) optical flares are a common phenomenon. In more 
luminous stars, flares are usually only detected through 
UV or X-ray observations ( Doyle et al. 1989b| ), although 
optical flares have been detected in young early K dwarfs 
like LQ Hya ( Montes et al. 1999| ). 

One would like to trace all the energetic processes in 
a flare to some common origin likely to be localized at 
some unresolvably fine scale in the magnetic field. The 


dicate that heating may propagate down to low photo- 
spheric levels, with densities higher than 10 16 cm -3 , al¬ 
though electrons with energies in the MeV range would 
be required to attain such depth. 

As regards the derivation of physical parameters such 
as electron density and temperature, various methods 


have been used. For example, Katsova (1990) published 
an analysis of flare Balmer decrements. The broadening 
and merging of higher Balmer lines, dominated by the 
Stark effect, have also been used to estimate electron den¬ 


sities in the chromosphere ( Donati-Falchi et al. 1985 ). 
Measurements of the broadening of the lower Balmer lines, 
which are less affected by the Stark effect, together with 
line shifts, provide information on the large scale mo¬ 
tions during flares. Steep decrements are evidence for 
electron densities between 10 8 and 10 12 cm -3 , while a 
shallow Balmer decrement indicates densities larger than 
10 13 cm~ 3 in the flaring M dwarf atmosphere. 
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Jevremovic et al. (1998) and Jevremovic (1999 ) devel¬ 
oped a procedure to fit the Balmer decrement based on 
the solution of the radiative transfer equation using the 
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esca pe probability comput ing technique of Drake (1980| ) 
and Drake & Ulrich (1980|) and the direct search method 
of |Torczon (1991| , |1992| ). 

In this paper we present the results of a large flare 
observed on AT Mic (Gliese 799B, V=10.25, B=11.83 
20' l 41 m 5U,-32°26'07", Equinox 2000 ) It is a visua l bi¬ 
nary star with a separation of 4.0 arcsec |Wilson (1978 ), ap¬ 
proximately 10.2 pc away, with both components of spec¬ 


tral type dM4.5e and subject to flaring (Joy & Wilson 
1949|)~ Kunkel (197C| ) determined the flare incidence at 
2.8 flar es per hour. Fr om its position in the (Mv, R-I) di¬ 
agram, Kunkel (1973|) deduced that AT Mic lies slightly 
above the main sequence, pointing t o its probable rnern- 
bership of the young disc population. Nelson et al. (1986 ) 
reported 29 U-band/B-band flares in 39.3 hours, only one 
of which showed microwave flaring at 5 GHz. They also 
reported two microwave bursts without any optical coun¬ 
terpart. Kundu et al. (1987 ) found that both components 
of the AT Mic system were active and variable at 6 and 20 
cm wavelengths. The peak emission levels were detected in 
the southern component (Gliese 799B). X-ray and ultravi¬ 
olet emission from this system has been observed several 
times over the last ten years or so ( |L insky et al. 19~82 ; 
Pallavicini et al. 199C). 

In §2 we present the current dataset plus a discus¬ 
sion of the data analysis while §3 gives the results. In §4, 
we analyse the evolution of the main physical parameters 
during the flare using the code developed by Jevremovic 
et al. ( 1998P and apply the procedure to a stellar and a 
solar flare for which parameters have been obtained using 
other techniques. The conclusions are given in §5. 


2. Observations and Data Reduction 

The observations were carried out at the South African 
Astronomical Observatory between 20-26 August 1985. 
Optical spectroscopy was obtained with the 1.9m 
telescope at SAAO equipped with the RGO Unit 
Spectrograph and the Reticon Photon Counting System 
(RPCS), providing spectral coverage from 3600 to 4600A 
at a reciprocal dispersion of 50A/mm and a wavelength 
resolution of AA ~ 2A (FWHM) at Hy. The spectral re¬ 
gion includes the entire Balrner series except Ha and H/3, 
as well as the Ca n H & K lines and the He I 4026A and 
He I 447lA lines. The exposure time was 60 sec. which 
with readout overheads gave a time resolution of 70 sec. 
During the run on AT Mic, sky conditions were generally 
good although seeing was variable. 

Spectra were reduced in the standard manner, con¬ 
sisting of background subtraction and flat-field correction 
using exposures of a tungsten lamp. Wavelength calibra¬ 
tion was made with Cu-Ar lamp spectra taken immedi¬ 
ately preceeding and following the stellar spectra. A first- 
order spline cubic fit to some 16 lines achieved a nomi¬ 
nal wavelength calibration accuracy of 0.12 A. An atmo¬ 
spheric extinction correction for SAAO was also applied to 
the data. Flux calibration, by reference to spectrophoto- 
rnetric standard stars (LTT7987 and LTT9239), was car- 



il! 


Fig. 1 . Normalised flare-only light curves for H6, Can K 
and the continuum (3600-3700A). The demarcation be¬ 
tween the impulsive and gradual phases, taking into ac¬ 
count the different line behaviour, is indicated by a vertical 
line. 


ried out using the NO AO packages IRAF.[] The observa¬ 
tions were made in spectrophotometric mode with a slit of 
400 /r width (~ 4 arcsec), adequate to encompass the en¬ 
tire image during the worst seeing conditions encountered 
during the run. A dekker was selected which provided an 
approximately square aperture. The programme star (AT 
Mic) and the spectrophotometric standards were observed 
under similar conditions. From a comparison of the cali¬ 
brations afforded by the two standards, we estimate the 
spectrophotometric accuracy to vary from ±0.12 magni¬ 
tudes (~ 12%) at 3600A to ± 0.03 magnitudes (~ 3%) 
at 4600A. The fluxes of individual emission lines will have 
additional errors associated with the photon noise level of 
the binned fluxes. 

3. AT Mic flare at 1:20 UT on 21 August 1985 

In this paper we report on an energetic flare, lasting more 
than one and a half hours, observed on 21 August 1985. 
The continuum emission, which initially peaked in the im¬ 
pulsive phase at 1:15 UT, had a second peak at 1:20 UT, 
roughly coincident with the primary peak in the Balrner 
lines. In addition, there is evidence for a second peak at 
1:30 UT in some of the Balrner lines (HA, H8) and also 
in He I 4026A. Fig. 0 shows a rapid continuum flux rise 
and decay, while the hydrogen Balrner lines peak later and 
decay more slowly. There is some indication that the low- 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper¬ 
ative agreement with the National Science Foundation. 
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order Balmer lines decay more slowly than the high-order 
lines, as has been noted by Butler (1991) in a flare on YZ 
CMi. The Can K line rises and decays even more slower 
than the hydrogen lines. Following Hawley & Pettersen 
(199l|J] we define the start of the gradual phase as the 
time when the continuum shows a turnover from fast to 
slow decay. The change from impulsive to gradual phase 
is indicated in Fig. |l] with a solid vertical line. 

Many emission lines other than the Balmer line series 
and Call H & K were seen during this flare, e.g. Fei and 
He I. Many of these were visible only during the impulsive 
phase and at the beginning of the gradual phase. Observed 
at higher spectral resolution these lines would provide im¬ 
portant information on the changing atmospheric struc¬ 
ture and hence flare heating. Fig. |j illustrates a time se¬ 
quence of the optical flare-only spectra obtained during 
the AT Mic flare. 

In order to form a mean quiescent spectrum we have 
co-added 75 quiescent spectra. The Balmer lines from Hy 
to H13 and the Can H & K lines are visible and in emis¬ 
sion in the quiescent spectrum. In the flare spectrum the 
Balmer lines increase their emission and there is a sub¬ 
stantial enhancement of the ultraviolet continuum. There 
is also some evidence for a filling in of the strong absorp¬ 
tion line Cal at 4227A observed during the impulsive and 
gradual phase of the flare. The He I 4026A and He I 447lA 
lines are also in emission during the flare. 

Emission line fluxes were measured as a function of 
time for the important contributors to the radiative losses 
in the optical, including the hydrogen Balmer lines from 
Hy to H9, the Call H & K lines and the Helium lines (see 
Table 1). To analyze the optical spectra, the mean quies¬ 
cent spectrum was first subtracted from the spectra taken 
during the flare. The emission line fluxes were then deter¬ 
mined by numerically integrating the flux under the line, 
using carefully chosen windows to set the local continuum 
level. Our low resolution spectroscopy cannot resolve the 
Can H and He line, but the blend is dominated by the 
behavior of He. We note the comparatively slow response 
of the Can K line compared to the Balmer lines. 

Bo th the Balmer lines and Can ('mission lines were 


i n solar flares is generally attributed to Stark broadening 
( Svestka 1972 ). 

In order to obtain an estimate of the broadening of 
the hydrogen lines during the flare, we have measured the 
full width at half maximum (FWHM) for Hy, H5, H8 and 
Can K. The FWHM for the Balmer lines were obtained 
using a Lorentzian profile, while a Gaussian profile was 
used for Can K. From Fig. |] we can conclude that the 
FWHM of the emission lines increases dramatically during 
the flare for the Balmer lines, although for Can K there 
is almost no change. 

It seems that line broadening continues well into the 
gradual phase of the flare, finishing around the same time 
as the continuum enhancement. As Doyle et al. (1988 ) 
pointed out, a possible explanation for the line wing en¬ 
hancement could be motions of several hundred kms -1 . 
If that were the case in the AT Mic flare, it would re¬ 
quire the amount of radiating material moving upwards 
and downwards to be approximately equal throughout the 
flare in order to produce the observed profiles. Although 
this could be the case it has not been observed on the Sun, 
where normally a pronounced redshift in Ha is usually ob¬ 


served during the impulsive phase (Zarro & Canfield 1989 
Wuelser fe Marti 1989). 


If hydrodynamic effects were contributing significantly 
in the chromospheric line forming regions, we would ex¬ 
pect Can K to show similar broadening to the Balmer 
lines, but this was not observed. Although the broadening 
of the Balmer series lines decreases almost exponentially 
during and after the impulsive phase, the central line in¬ 
tensity of HJ (Fig. []]) shows two peaks, indicating that 
there are two energy releases. The first burst induced line 
broadening while the second burst did not, suggesting that 
this flare may have derived from an arcade of loops. 

The excess broadening of the Balmer lines could be 
interpreted as turbulence as we did not observe any clear 
asymmetry at any time, in any line, during the flare. 

Our conclusions in the last two sections, regarding 
the evolution and broadening of the emission lines dur¬ 
ing the flare on AT Mic, are similar to those of Hawley & 


above hhc quiescent values for more than 2 hours, while 
the continuum enhancement lasted around one hour. 


3.1. Line Broadening 

Near flare maximum, the wings of Ha, H<5, H8 and even 
H9 are very broad, around 25A at the base, as is shown 
in Fig. [| As a result, the Balmer lines start to merge at 
about H13; meanwhile Can K shows very little broad¬ 
ening. piampapa (1983 ) and Worden (1983| ) pointed out 
that during flares the lower Balmer lines such as Ha and 
H/3 were unbroadened, while the higher Balmer series lines 
show significant broadening, 20A or more for some lines. 
They also reported that the Balmer lines (and Hei) in 
a very large flare were broadened and asymmetric to the 
red, while Can is unbroadened. Balmer line broadening 


Pettersen (1991) for a flare on AD Leo. 


3.2. The He I 4026A and He I 4471 A Triplet Lines 

Helium lines are produced under generally higher exci¬ 
tation conditions than other lines in the visible and the 
lines are generally optically thin, which makes interpreta¬ 
tion easier. The helium lines are also useful as indicators 
of the presence of electric fields in the flaring plasma due 
to their forbidden components. However, even for the Sun, 
exact knowledge of the background spectrum is required 
for the determination of the contribution of the forbidden 
lines. The He I lines have been reported in solar white light 
flares (bites et al. 1986) and i n stellar flares (Moffett fe 
Bopp 1971 ; Doyle et al. 1988|) . Lites et al. (1986 ) pointed 
out that any mechanism which produces a dense plasma at 
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Fig. 2. The time sequence of optical flare-only spectra of AT Mic during the pre-flare phase, impulsive phase, maximum 
continuum, maximum flux in Hd, minimum between the two maxima, second maximum flux in H<5, the beginning of the 
gradual phase and finally the gradual phase. The mean quiescent spectrum is shown as a dotted line for comparison. 


Table 1. The total energy radiated in the Balmer, Calcium and Helium lines during the flare on AT Mic. Units are 
10 30 erg 


UT 

h 7 

H<5 

H8 

H9 

H10 

Can I< 

Can H+He 

Cal 

He I 4026 A 

Hei 4471 A 

1:13 

1.61 

1.42 

.74 

.46 

.24 

.04 

.36 

.04 

.52 

.90 

1:15 

4.62 

4.37 

2.50 

2.10 

.87 

.25 

3.76 

.01 

1.34 

1.43 

1:20 

11.07 

10.69 

7.54 

5.28 

2.83 

2.28 

9.86 

.36 

.67 

1.24 

1:25 

9.03 

8.25 

6.12 

5.24 

3.23 

2.42 

8.28 

.14 

.56 

.94 

1:28 

8.12 

8.52 

6.32 

4.91 

2.84 

2.60 

8.10 

.41 

.61 

.90 

1:30 

6.30 

5.22 

4.05 

3.42 

2.08 

1.63 

5.94 

.34 

.48 

.48 

2:00 

1.81 

1.55 

.97 

.85 

.65 

.51 

1.85 

.11 

.13 

.39 


relatively high temperatures (2 10 4 — 10 5 K) will produce 
emission in neutral helium. 


The He I 4026A line is a diffuse triplet between levels 
2 3 P and 5 3 D. This line is usually not detected except in 
the strongest events. It was shown (e.g. Vidal 1964) that 
the higher members of this triplet are quite sensitive to 
density changes. The He I 4026A line is near to a strong 
Mnl absorption feature at 4030A (which during intense 
act ivity can go into emissio n due to pumping from Mgn 
k, ( Doyle et al. 1992 , 200 1| )) thus making detection diffi¬ 
cult. We observed the He I 4026A line and Hei 4471 A line 


in emission during the impulsive phase and well into the 
gradual phase of the flare on AT Mic (Fig. |j). Two peaks 
are seen in the time profile of Hei 4026A, the first at 1:20 
UT and the second at 1:30 UT. These are simultaneous 
with the maxima in the Balmer lines. In Hei 447lA only 
a single peak at 1:20 UT is visible. 


The FWHM of the two Hei lines are 5.0 A and 2. 8A 
respectively. Gieske fe Griem (1969 ) give a procedure 
for computing the line profile of Hei 4026A including 
Stark and Doppler broadening. With their model for 
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Fig. 3. Comparison of the mean quiescent and flare-only line profile for the Balmer series lines, Can H + He and Can 
K. 


T e = 20 000 K, we derive an upper limit for the electron 
density of 10 16 cm -3 for the He I emitting region. 

3.3. Continuum Analysis 

We have measured the fluxes (Table 2) for three differ¬ 
ent continuum regions: one in the near ultraviolet, 3600- 
3700A and two in the blue, 4130-4300A (with the ex¬ 
clusion of the Cal line region) and 4435-4545A. Fig. 
shows the variation in the ratios of the continuum fluxes 
during the impulsive and gradual phases of the flare. We 
note that the ratios of the near ultraviolet band to the 
two longer wavelength bands do not reach their maxi¬ 
mum blueness until the end of the impulsive phase and 
the beggining ot the gradual phase (indicated by the verti¬ 
cal line), whereas the ratio of the intermediate wavelength 
band to the longer wavelength band reaches it maximum 
blueness at the beginning of the impulsive phase. As all 
three continuum bands reach their maximum brightness 
at the start of the impulsive phase, this behaviour suggest 


that a simple heating of a static, confined, kernel with ra- 
diation follo wing a black-body law, is not appropiate (see 
Butler 1996). 

hJopp & IViottett (19Y3) suggested that the continuum 
enhancement ends when the fast decay (what we term 
the impulsive phase) has ended. In fact, this generally ac¬ 
cepted notion, that the white-light continuum is present 
only during the impulsive phase is certainly not the case 
for the AT Mic flare nor for th e flare on AD Leo reported 
by Hawley fe Pettersen (1991 ). 

Spectroscopic observations of M dwarf flares reveal the 
shape of the flare continuum and allow us to explore the 
possible continuous emission mechanisms. It has been sug¬ 
gested (e.g., Cram & Woods 1982; Houdebine 1992) that 
the optical flare continuum is composed of bound-free and 
free-free emission from a two-component thermal model 
characterized b y electron temperatures of T ~ 10 5 K and 
T ~ 2 10 4 K. Cram fe Woods (1982 ) showed the spec¬ 
tral energy distributions of some of the suggested radia¬ 
tive processes that may give rise to the flare continuum. 
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Fig. 4. FWHM through the flare-only spectrum on AT 
Mic for Hy, H5, H8 and Can K lines. For the Balmer lines 
we have used a Lorentzian profile and Gaussian profile for 
the Ca II K line. In the top left we plot the Voigt profile 
fitting (dotted) for Hy. 




4020 4025 4030 4035 4450 4460 4470 4480 4490 

Wavelength (A) Wavelength (A) 

Fig. 5. Left panel: Average over 4 flare-only spectra of 
He I 4026A around flare maximum. Right panel: Hei 
447lA average over 4 flare-only spectra around flare max¬ 
imum. 


The shape of the free-free curve at T ~ 10 5 K obtained by 
Giampapa (1983) produces reasonable fits to the observed 
flare continuum. Kunkel (1970) has argued that hydrogen 
free-bound emission is the principal component for the ob¬ 
served continuum emission in M dwarf flare ev ents. The 
curve for this process given by Giampapa (1983 ) only fits 
the red part of the continuum, which would imply that 
additional continuum processes, other than bound-free, 


are taking place in higher temperature layers. Hawley & 
Fisher ( 1992|) fitted the observed flare continuum, but the 


Fig. 6. Continuum flux ratios through the flare on AT 
Mic. It is clear from the plot that the continuum of the 
flare has a sustantial blue component during both the im¬ 
pulsive phase and gradual phase and that the maximum 
blueness in the near UV occurs close to the end of the 
impulsive phase. 

Table 2. Continuum energy radiated during quiescence 
and during flare on AT Mic. Units are 10 32 erg 


UT 

3600 - 3700 A 

4150 - 4300 A 

4435 - 4545 A 

Quiescence 

.03 

.05 

.12 

1:13 

.04 

.05 

.11 

1:15 

.98 

.68 

.74 

1:20 

.63 

.31 

.35 

1:25 

.40 

.17 

.22 

1:28 

.37 

.16 

.21 

1:36 

.23 

.11 

.16 

2:00 

.06 

.06 

.12 


results were much bluer than those computed from the 
models. They concluded that the best fit to their observa¬ 
tions was a blackbody spectrum with T ~ 8500 — 9500 K. 
[Hawley fc Fisher (1992 ) proposed that the flare contin¬ 
uum is formed by photospheric reprocessing of intense 
UV/EUV line emission from the upper chromosphere. 

Game & Vial (1985) proposed a model for solar white- 
light flares in order to explain the observed blue contin¬ 
uum. Using the study of radiative processes in solar flares 
( Dame fc Cram 1983| ) they show how existing semiem- 
pirical models can be modified to predict a visible spec¬ 
trum more similar to the WLFs emission observed. The 
hot temperature pl ateau (6-9 10 4 K) model proposed by 
Dame fe Vial (1985 ) is not valid for all WLFs but provides 
an explanation for the origin of an unidentified source of 
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opacity distinct from the Balmer and Paschen continuum. 
Houdebine et al. (1996| ) suggested that this blackbody 
type emission is due to the enhanced electron density in 
the lower chrosmophere. 

However, the blue continuum enhancement could also 
be due to merging of the higher members of the Balmer 
lines particularly in the vicinity of the Balmer jump, 
which indicates that these lines contribute a non-negligible 
amount to the near-ultraviolet continuum (Donati-Falchi 
Zarro fe Zirin 1985| ). One en re- 


et al. .985 


ported in which the continuum was found to be red in c olor 
during the first stages of the flare (|Doyle et al. 1989a ). 


the best solutions for the Balmer decrement, the procedure 
allows us to calculate the effective thickness of the slabs 
of hydrogen plasma with which we are able to determine 
the surface area of the emitting plasma from the calcu¬ 
lated emission measure. Although the temperature of the 
underlying source could be a free parameter in our code, 
we have fixed a lower limit at 5 500 K for the Sun and 
2 500 K for YZ CMi and AT Mic, for numerical stability 
reasons. 

4.2. Comparison with other Codes 


4. Derivation of Flare Parameters 

We normalized the flare-only fluxes for Hd, H8, H9 and 
H10 to Hy at several times during the impulsive and grad¬ 
ual phases. We do not include Hll and higher Balmer lines 
due to the difficulty in assigning the local continuum level. 
The Balmer decrement became quite shallow and stayed 
nearly constant during the impulsive phase of the flare 
and well into the gradual phase. After this, the decrement 
steepens but does not become as steep as that during the 
quiescent phase. 

4.1. Our Procedure 

Jevremovic et al. (1998) developed a procedure to fit the 
Balmer decrements. It is based on the solution of the 
radiative transfer equation, using the escape probability 
technique of Drake (1980|) , and Drake fe Ulrich (1980p . 
The Drake fe Ulrich (1980| ) method assumes: (i) that the 
escape probability for a slab of optical thickness is av¬ 
eraged over all directions with a uniform distribution of 
emitters, (ii) complete redistribution of frequencies and 
(iii) that the line profile is represented by a Doppler core 
and p ower law wings of slope — 5/2 to emulate the Stark 
effect. Jevremovic et al. (1998 ) adopted a simplified pic¬ 
ture of the flaring plasma as a slab of hydrogen with an 
underlying thermal source of radiation and which causes 
photoionization. This source, which can be considered to 
be the photosphere, can be exposed to additional heat¬ 
ing during the flare. It assumes that the Balmer lines are 
formed in the upper chromosphere. We note that, this sim¬ 
plified picture leaves out some important processes (i.e. 
cooling from MgII, Can and other metals). However, as 
seen from Table 1 the AT Mic flare is predominantly seen 
in Balmer lines and we believe that this picture may be ap¬ 
plicable. The procedure minimizes the difference between 
the observed and calculated Balmer decrement using a 
multi-directional search algorithm (Torczon 1991, 1992Q . 
This allows us to find the best possible solution for the 
Balmer decrement in four parameter space where the pa¬ 
rameters are: electron temperature, electron density, op¬ 
tical depth in the Lya line and the temperature of the 
underlying source. Details of the procedure and the multi¬ 
search algorithm and com parison with the classic al sim¬ 
plex algorithm are given in Jevremovic et al. (1998 ). From 


In order to check the |Jevremovic et al. (1998| ) procedure 
we have applied it to two flares for which physical param¬ 
eters have been independently derived, one stellar and one 
solar. 


4.2.1. Comparison with the Gas-Dynamic Model for 
the YZ CMi flare 

We have applied our procedure to a flare on YZ CMi ob¬ 
served during a coordinated run involving telescopes at the 
South African Astronomical Observatory, and the ESA 
X-ray satellite, EXOSAT, on 4 March 1985. The photo¬ 
metric observations consisted of UBVRI photometry with 


a 0.75-m telescope. Details have been published by poyle 
|et al. 


. The spectra of the flare on YZ CMi, over the 
wavelength range 3600-4400A were analysed by Katsova 
|et al. (1991 ) using a gas-dynamic model. Based on the 
gas-dynami c model and using the data se t for the flare 
observ ed by poyle et al. (198gj ) on YZ CMi, Katsova et al 


(1991 ) derived the main physical parameters at flare max¬ 


imum. They matched quite well the observed data with 
the theoretical decrement for a temperature T = 10 4 K, 
electron density 10 14 cm -3 and optical depth of the layer 
at the Lya line center T£ ya = 210 6 . They also derived 
the emitting source area of > 5xl0 17 cm 2 using a descrip¬ 
tion of the optical continuum as a Planck function for a 
temperature of 10 4 K. 

Our values for the derived flare p arameters are reason- 
ably close to those values obtained by Katsova et al. (1991 ) 
(see Table 3), although our values for the electron tem¬ 
perature and Lya optical depth at line center are slightly 
larger and smaller respectively. We obtained an effective 
thickness of the slab of hydrogen plasma of around 2 000 
km and an emissivity of 3.25 erg cm 


^ 3 s - 4 


4.2.2. Comparison with a Stark Broadening Model for 
a Solar flare 

As a second check we have also applied our procedure to a 
solar flare observed on 6 March 1993. Details of the flare 


have been published by |johns-Krull et al. (1997| ). They 
fitted the observations with models that include Stark 
broadening by electrons, assuming that the Balmer lines 
form in a homogeneous flaring region overlying the quiet 
photosphere that has n e « (2 — 5) 10 13 cm -3 . They mod- 
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Table 3. Comparison of physical parameters for flares on YZ CMi and on the Sun obtained with the different codes. 


Star 

Code applied 

Reference 

N e (cm" 3 ) 

Te( K) 

r{Lya) 

Area(cm 2 ) 

YZ CMi 

Our Code 

1 

5 10 13 

13 000 

310® 

no 18 


Gas Dynamic 

2 

10 14 

10 000 

210® 

5 10 17 

The Sun 

Our Code 

1 

2 10 14 

18 000 

7 10 4 

710 19 


Stark Broadening 

3 

(2-3) 10 13 

14 000 

n/a 

5 10 19 


References'. 1: This work, 2: Katsova et al. (1991), 3: Johns-Krull et al. (1997). 


elled the line widths following the g eneral procedure out - 
lined by Redman fc Sucmoto (1954 ) and [Svestka (1963|) , 
which attempts to model the flare emission profile by tak¬ 
ing into account the effects of self-absorption, Doppler 
broadening and Stark broadening, and comparing the re- 
sulting profiles with the observed flare-quiet Sun profiles. 
Johns-Krull et al. (1997) found that, over the course of 
the first 10 flare spectra, the electron density varied fairly 
randomly in the range (2 — 3) 10 13 cm -3 . Meanwhile, the 
values for the electron temperature were around 14 000 K. 
Although they did not calculate the flare area they pro¬ 
vided an Ha image taken at BBSO around 10 minutes 
after the flare maximum and from this image one can in¬ 
fer a solar flare area of about 5 10 19 cm 2 . 

Using our procedure, an d the solar flare data obtained 
by Johns-Krull et al. (1997| ), we obtain the parameters at 
solar flare maximum as shown in Table 3. As in the case of 
the YZ CMi flare, our value s are close to those obtained 
by Johns-Krull et al. (1997 ) although our value for the 
electron density is almost one order of magnitude larger. 
This result for the electron density could be due to the fact 
that, in our fit of the solar flare, we have not fitted the 
lower Baimer lines, Ha and H/3 which are formed higher in 


the atmosphere, and were observed by Johns-Krull et al. 
(199?j)] We have also discarded H7 in our fitting, in order 
to compare the results of both solar and stellar flares based 
on the same Balmer lines. The flare area is in very good 
agreement with that derived from the observations. We 
obtained an effective thickness for the hydrogen slab of 
almost 100 km and an emissivity of 145 erg cm -3 s _1 . 


4.3. Modelling the AT Mic flare 


From the above two cases we concluded that the procedure 
of Jevremovic et al. (199S| ) gives reliable results and we 
proceed to apply it to the AT Mic flare. In Fig. [7] we have 
plotted the data and three Balmer decrement fits during 
the flare and one during quiescence. The evolution of the 
main parameters is given in Fig. || 

Analysing the Balmer decrement during the AT Mic 
flare, we observed that, at the beginning of the impulsive 
phase, the Balmer decrement has a steep slope which be¬ 
come shallower as the flare progresses through its impul¬ 
sive phase. The Balmer decrement keeps almost the same 



Fig. 7. The observed Balmer decrements (squares) and 
the optimum computed fits (solid line) during quiescence 
and at different times during the flare. 


slope until well into the gradual phase, after which the 
slope steepens reaching similar values to that in the qui¬ 
escent phase. In the two spectra prior to the first H(5 max¬ 
imum at 1:20 and in the three spectra after the second H(5 
maximum at 1:28, the H(5 flux is larger than the Hy flux. 
Due to their similarity, we have compared the observed 
Balmer decrement behaviour for the flare on AD Leo re¬ 
ported by Hawley & Pettersen (1991) and the observed 
Balmer decrement for our flare on AT Mic. During the 
quiescent phase, the slope for AT Mic is steeper than for 
AD Leo, while, during the impulsive and gradual phases, 
the slopes for AT Mic are shallower than thoses for AD 
Leo. This means that the variations of the physical pa¬ 
rameters during the AT Mic flare are larger than during 
the AD Leo flare. 

From Fig. H one can see that there are three distinct 
peaks in the electron temperature, of which two coincide 
with peaks in HA We interpret these two peaks as the re¬ 
sult of two sudden separate releases of energy (which pre¬ 
sumably correspond to the two continuum peaks shown 
in Fig. (lj) , while the third peak in the temperature is con- 
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Fig. 8. Parameters derived from the fit to the Balmer decrement using the direct search method. Evolution of the 
electron temperature and temperature of background source are shown in the top left panel, electron density in bottom 
left panel , effective plasma thickness in top right panel and goodness of fit in bottom right panel. 


nected to additional, more gradual heating. This is also 
seen in the behaviour of the temperature of the back¬ 
ground source and in Fig. £]. Note that peaks in the slab 
thickness coincide with the maximum temperature, which 
can be interpreted by a change in the degree of ionization 
with temperature. An increase in the degree of ionization 
leads to a thicker slab for the same optical depth. Also, 
the thickness of the slab, coupled with the low emissivity 
of the plasma, explains the relatively low increase in the 
continuunr/line flux during the third peak (see Fig. [I]). For 
example, the modeled emissivity in Hy for the first peak 
was around 150 erg cm -3 s -3 while during the third event 
the emissivity was around 10 erg cm -3 s -3 . Similar vol¬ 
umes of plasma for both peaks lead to much lower increase 
in the flux during the third event (contrary to what one 


would expect taking into account only thermal radiation). 
Also, one has to keep in mind that our background source 
is included only as the source of photoionization which 
controls ionization balance. Inclusion of other sources of 
ionization (i.e. focused beam of electrons/protons) would 
probably change the results of modeling, but that is be¬ 
yond the scope of this paper. 

We found little change in the electron density during 

the flare(s) and it stays between 10 13 -10 14 cm -3 . These 

values are one or two orders of magnitude below the up¬ 
per limits deduced from either Stark broadening or the 
Inglis-Teller relation. This confirms the reliability in our 
method, based in the simultaneous analysis of a few lines, 
rather than some earlier estimates of electron density (i.e. 

, based on a single line analysis. 


Houdebine 1992 
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From the best fit solutions for the Balmer decrement, 
we have calculate the effective thickness of the slab of hy¬ 
drogen plasma from the beginning of the impulsive phase 
until the beginning of the gradual phase as shown in the 
top right panel of Fig.||. Using these results we are able 
to determine the surface area of the emitting plasma from 
the calculated emission measure and the measured fluxes 
in Hq. From these results we conclude that the area of 
the flare reach three maximum values, which take place 
after each maxima in the electron temperature. The val¬ 
ues we obtain are around few times 10 19 cm 2 for all of 
them. Taking into account that the radius of AT Mic is 
0.32 Rq, this suggests that the flares occupy around 0.6 
percent of the stellar surface area. This is comparable with 


the siz es of t he largest solar flares (|Tandberg-Hanssen fc 
Emslic 1988 ). 


5. Conclusions & Discussion 

We observed a major flare on AT Mic in the 3600-4500A 
wavelength region which allowed us to compute a detailed 
flare energy budget in the optical spectral region. The 
physical parameters obtained are consistent with previ¬ 
ously derived values for stellar flares. We have compared 
the results from our method for two flares previously anal¬ 
ysed, namely, one on YZ CMi and one on the Sun, obtain¬ 
ing reasonable agreement. 

We have obtained the first detailed trace of physical 
parameters during a stellar flare. We have also shown that 
our method could provide a suitable model for the study of 
flares on stars of different spectral type (G-M). The total 
radiated energy released in the region 3600-4500A during 
the flare was ~ 4 1 0 33 erg and A U ~ 4 mag. The energy 
radiated by this flare is comparable with the strongest 
solar ones but is of medium strength compared to flares 
on other dMe stars. 

Further stellar observations with high time and high 
spectral resolution should be obtained. Also, further so¬ 
lar observations should be obtained in order to check the 
method with other examples, particularly, the flare area 
and thickness. 
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